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Dielectric  Resonators  for  Use  in  Microwave 
and  Millimeter-Wave  GaAs  Circuits 


1 .  Background 

Dielectric  resonators  are  of  considerable  interest  for  use  in  higher 
microwave  and  millimeter-wave  frequency  stabilization  circuits  and  filtering 
structures.  As  frequency  increases  dissipation  losses  in  conducting  materials 
make  microstrip-like  circuit  elements  too  lossy  (low  Q)  and  metal  waveguide 
structures  are  expensive  to  machine.  Also,  as  frequency  increases,  lumped 
element  structures  become  difficult  to  realize  and  are  also  low  in  Q. 
Although  dielectric  resonators  have  been  of  interest  for  some  time  [1,2], 
little  effort  has  been  placed  on  investigating  dielectric  resonators  which  are 
suitable  for  processing  in  monolithic  microwave  and  millimeter-wave  circuits. 
It  should  be  noted  that  recently,  some  work  has  been  described  [3]  for 
coupling  monolithic  circuits  to  non-mono lithic  dielectric  resonators.  Since 
GaAs  is  a  substrate  material  receiving  considerable  attention  for  processing 
monolithic  microwave  and  millimeter-wave  circuits,  this  report  describes  the 
preliminary  phase  of  an  investigation  into  dielectric  resonators  which  might 
be  capable  of  being  processed  in  or  on  GaAs  to  interface  with  surrounding  mono¬ 
lithic  circuits. 
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2. 


Theoretical  Formulation 


The  initial  dielectric  resonator  structure  considered  for  use  in  GaAs 
circuits  is  depicted  in  Fig.  1.  Figure  la  shows  a  perspective  sketch  of  such 
a  resonator,  while  Fig.  lb  shows  a  cross  section  through  the  resonator  and 
defines  various  parameters  associaed  with  the  geometry.  It  should  be  noted 
that  H  (height  of  the  groove  bottom  above  the  lower  ground  plane)  is  taken  to 
variable  down  to  H  equal  to  zero,  where  the  dielectric  resonator  is  sitting  on 
the  lower  ground  plane.  For  a  given  thickness  of  GaAs  substrate  Tg ,  T  and  H 
are  related  by 

Tg  =  T  +  H  (1) 
The  initial  concept  behind  this  kind  of  resonator  is  to  determine  if  reason¬ 
able  geometries  for  use  at  approximately  95  GHz  can  be  found  where  the  groove 
could  be  etched  into  GaAs  by  "dry  etching"  methods  (e.g.  reactive  ion  etch¬ 
ing)  . 

2.1  No  Fringing  Fields 

In  order  to  estimate  the  geometries  feasible  for  this  structure  at 
95  GHz  an  approximate  analysis  was  developed.  This  analysis  employs  an 
approach  similar  to  that  used  by  M.  Dydyk  [4]  for  dielectric  resonator 
cylinders  placed  upon  some  substrate  material  (usually  different  from  that  of 
the  resonator  material)  and  coupled  to  a  microstrip  line.  The  initial  treat¬ 
ment  employed  for  the  work  in  this  report  approximated  the  resonator  of  the 
type  shown  in  Fig.  1  by  assuming  that  the  groove  is  very  wide  and  thereby 
enabled  the  Fig.  1  resonator  to  be  analyzed  by  the  treatment  described  in  [4 
and  5].  This  solution  is  based  upon  a  transverse  resonance  method  described 
in  [6].  This  approach  assumes  a  magnetic  wall  at  the  cylinder  boundary  and 
thereby  neglects  fringing  fields.  The  effects  of  these  fringing  fields  will 
be  accounted  for  in  a  refinement  to  the  analysis  as  shown  later. 
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Plane 


g.  la  —  Perspective  sketch  of  groove-type  dielectric  resonator  in  GaAs. 
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Fig.  lb  —  Cross  section  of  groove-type  dielectric  resonator. 
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For  the  no  fringing  case  the  expression  which  interrelates  the  geometry 
and  material  parameters  for  this  approximate  representation  is 

T  _  1 -  -1  ,  YoaR  d,H., 

2R  "  2y.R  ^tan  ^  y.R  Coth  YoaR(R)] 


+  tan_1  [  C°th  YosR(l)1} 


In  (2), 

v  =  7(^) 2£  .  ( 


oa 


=  7(^1 ~2  -  (p 
=  J(~^>2  - 


(2) 

(3) 

(4) 

'OS  *  '  R  '  VC'  "s 

Here,  c  is  the  velocity  of  light  in  free  space,  u>  is  the  angular  frequency  of 
operation,  is  the  relative  permittivity  of  the  cylinder  dielectric  mate¬ 
rial,  £g  is  the  relative  permittivity  of  the  substrate  and  surrounding  mate¬ 
rial  on  which  the  cylinder  sits.  For  the  special  case  of  the  groove-type 
resonator  in  GaAs,  =  £g  =  13.  The  remaining  parameters  in  equations  (2) 
through  (5)  are  defined  in  Fig.  lb.  Equations  (2)  through  (5)  were  programmed 
in  general  form  for  digital  computer  usage.  This  program  is  listed  and 
described  in  Appendix  A  of  this  report  (see  program  DRES). 

2.2  With  Fringing  Fields 

As  described  by  Dydyk  [5]>  a  refinement  can  be  made  in  the  use  of 
equation  (2)  to  account  for  fringing  fields.  This  is  accomplished  by  solving 
(6) 


J  (P  ,R)  K  (a  R) 

f*  R  (  TTFTT  1  =  -«  R  [  A„C°0,  ) 


ci“  1  JxOciR) 


co  1  Kl(acoR) 


(6) 
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for  p^R  and  replacing  the  ratio  (2.405/R)  in  equations  (3)  through  (5)  by  the 
newly  found  p^  value.  This  process  accounts  for  the  decay  of  fields  external 
to  the  resonator.  In  equation  (6),  Jq  is  the  zeroth  order  Bessel  function  of 
the  first  kind;  J,  is  the  first  order  Bessel  function  of  the  first  kind;  K  is 
the  zeroth  order  modified  Bessel  function  of  the  second  kind;  and  is  the 
first  order  modified  Bessel  function  of  the  second  kind.  Also  in  equation  (6) 


a 

co 


2 


and 


Solving  (6)  through  (8)  to  eliminate  p  yields 


(7) 

(8) 


J  (0  -R) 

N  ~  r-  -i 


=  -Q  [ 


Kq(Q) 

K^Q) 


(9) 


(10) 


where  Q  =  J(^)2  (e^-D-O^R)2 
The  solution  of  equations  (9)  and  (10)  were  programmed  for  digital  computer 
usage.  This  program  is  listed  and  described  in  Appendix  B  of  this  report  (see 
program  BETACI). 

The  computer  programs  listed  and  described  in  Appendices  A  and  B  were 
tested  for  validity  by  analyzing  structures  employed  by  M.  Dydyk  [4) .  In  this 
previous  work  the  solutions  were  compared  to  experimental  results  and  were 
found  to  agree  to  within  less  than  two  percent  error  on  resonant  frequency. 
Figure  2  shows  plots  of  resonator  geometry  paramters  versus  frequency  (free 
space  wavelength)  for  a  fixed  substrate  thickness  to  radius  ratio  of  H/R  =  0.2 
with  variation,  in  the  height  of  the  upper  ground  plane  above  the  resonator 
(D/R).  Figure  3  shows  similar  data  for  the  case  (D/R)  -*  -with  H/R  varying. 
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(t/2R) 


Fig.  2  —  Simulation  of  design  geometries  for  circular  disk  dielectric 
resonators  on  a  dielectric  substrate  for  variations  in  upper  ground 
plane  spacing. 


(t  /  2  R ) 


Fig.  3  —  Simulation  of  design  geometries  for  circular  disk  dielectric 
resonators  on  a  dielectric  substrate  for  variations  in  substrate 
thickness. 
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For  both  Figs.  2  and  3  £  =  37.0  and  e  =2.1.  The  correlation  between 
Dydyk's  results  and  the  results  evaluated  with  the  programs  used  here  is 
extremely  close.  This  close  agreement  with  Dydyk's  computational  and  experi¬ 
mental  results  has  provided  verification  of  the  programmed  analyses  for  use  on 
other  configurations  within  the  scope  of  the  approximations  assumed  in 
treating  the  groove  and  coupling  mechanism  for  the  resonator. 

3 .  Design  Data 

The  analyses  described  in  section  2,  together  with  implementations  shown 
in  Appendices  A  and  B,  were  applied  to  the  case  where  £f  =  £g  =  13.  This  case 
is  representative  of  a  resonator  and  substrate  composed  completely  of  semi- 
insulating  GaAs .  Figures  4  and  5  display  plots  of  geometry  versus  frequency 
for  different  ranges  of  resonator  thickness  T  with  the  upper  ground  plane 
sufficiently  far  removed  (D/R  -*•  °°)  so  as  to  have  no  effect.  It  was  felt  that 
from  a  practical  standpoint  this  is  the  most  useful  case. 

From  this  data  it  was  intended  to  gather  two  sets  of  information.  First, 
design  data  was  to  be  generated  for  frequency  scaled  experiments  to  investi¬ 
gate  coupling  mechanism  (study  the  effect  of  various  groove  widths)  empiri¬ 
cally,  since  this  effect  is  not  treated  directly  by  the  analyses.  Scaled 
experiments  are  also  necessary  to  determine  the  separator  necessary  between 
the  resonator  and  microstrip  line  from  which  energy  is  to  be  coupled.  By 
scaling  to  lower  frequencies  (microwave)  more  accurate  measurements  could  be 
performed  with  readily-available  automated  network  analyzers.  Secondly, 
Figures  4  and  5  were  to  be  used  to  generate  designs  for  millimer-wave  reso¬ 
nators  to  determine  the  feasibility  for  geometries  in  GaAs  at  these  higher 


frequencies . 


(27TR/XJ 


'  0.1  0.2  0.3  0.4  0.5  0.6 

(t/2R) 


Fig.  5  —  Simulation  of  design  geometries  for  circular  disk  dielectric  (GaAs) 
resonator  on  dielectric  (GaAs)  substrate  for  this  substrates. 
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The  data  shown  in  the  following  tables  was  generated  to  be  consistent 
with  several  constraints.  The  data  for  frequency  scaled  work  (microwave 
frequencies)  was  generated  for  resonances  in  the  500  MHz  to  2  GHz  range.  For 
simplification  in  experiments  it  was  considered  desirable  to  work  on  TiO^- 
loaded  Stycast  (purchased  from  Emerson  &  Cumings,  Inc.)  sheets  having  K  =  13 
and  thicknesses  of  0.25  inch  and  0.50  inch.  This  would  enable:  (a)  resonator 
structures  to  be  fabricated  by  simple  machine  processes;  (b)  use  of  n-type 
connectors  for  launching  into  coupling  microstrip  lines,  and  (c)  use  of  metal 
(e.g.  aluminum)  foil  sheets  to  be  used  for  simulating  the  relatively  wide 
microstrip  lines  needed  for  use  at  these  frequencies.  Further  discussion  of 
these  experiments  will  be  provided  in  section  4  of  this  report. 

For  T  +  H  ^  0.250  inch,  D  ■+  ® 

Table  I.  (R  =  1.25  inch) 


(T/2R) 

(H/R) 

(2JtR/\  ) 

0 

T(IN) 

H(IN) 

f  (GH  ) 

0  z 

0.100 

0 

3.3 

0.250 

0. 

5.0 

0.075 

0.05 

3.4 

0.188 

0.062 

5.1 

0.050 

0.10 

3.7 

0.125 

0.125 

5.6 

Table  II. 

(R  =  0.625  inch) 

(T/2R) 

(H/R) 

(271R/X  ) 

0 

T(IN) 

H(IN) 

f  (GH  ) 

0  z 

0.2 

0. 

1.87 

0.250 

0. 

5.6 

0.175 

0.05 

1.88 

0.219 

0.031 

5.7 

0.150 

0.10 

1.92 

0.188 

0.062 

5.8 

0.100 

0.20 

2.2 

0.125 

0.125 

6.6 

0.050 

0.30 

2.8 

0.062 

0.188 

8.4 

l 
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Table  III. 

(R  =  0.3125 

inch) 

(T/2R) 

(H/R) 

(271R/A  ) 

0 

T(IN) 

H(  IN) 

f  (GH  ) 

0  z 

0.400 

0. 

1.17 

0.250 

0. 

7.0 

0.375 

0.05 

1.17 

0.234 

0.0156 

7.0 

0.350 

0.  10 

1.17 

0.219 

0.0312 

7.0 

0.300 

0.20 

1.20 

0.188 

0.0625 

7.2 

0.200 

0.40 

1.37 

0.125 

0.125 

8.2 

0.100 

0.60 

1.9 

0.0625 

0.188 

11.4 

0.050 

0.70 

2.5 

0.0312 

0.219 

15.0 

Table  IV.  (R  =  2.50  inch) 

(T/2R) 

(H/R) 

(2rtR/A.  )  T(IN) 

0 

H(IN) 

f  (GH  ) 
o  z 

0.05 

0. 

4.2  0.25 

0. 

3.2 

Table  IV. 

(R  =  5.0  inch) 

(T/2R) 

(H/R) 

(2nR /Kq) 

T(IN) 

H(IN) 

f  (GH  ) 

0  z 

0.025 

0. 

7.3 

0.25 

0. 

2.7 
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For  T  +  H  £  0.125  inch,  D/R  ■*  * 


Table  V. 

(R  =  0.625  inch) 

(T/2R) 

(H/R) 

(2nR/A  ) 

0 

T(IN) 

H(IN) 

f  (GH  ) 

0  2 

0.1 

0. 

3.3 

0.125 

0. 

9.9 

Table  VI. 

(R  =  0.3125 

inch) 

(T/2R) 

(H/R) 

(271R/A  ) 

0 

T(IN) 

H(IN) 

f  (GH  ) 

0  z 

0.200 

0. 

1.88 

0.125 

0. 

11.3 

0.175 

0.05 

2.05 

0.109 

0.016 

12.3 

0. 150 

0.10 

2.29 

0.094 

0.031 

13.8 

0.100 

0.20 

3.00 

0.032 

0.063 

18.0 

Table  VII. 

(R  =  0.15625  inch) 

(T/2R) 

(H/R) 

(27XR/A  ) 
o 

T(IN) 

H(IN) 

f  (CH  ) 

0  z 

0.400 

0. 

1.17 

0.125 

0. 

14.1 

0.375 

0.05 

1.17 

0.117 

0.008 

14.1 

0.350 

0.10 

1.17 

0.109 

0.016 

14. 1 

0.300 

0.20 

1.20 

0.094 

0.031 

14.4 

0.250 

0.30 

1.27 

0.078 

0.047 

15.3 

0.200 

0.40 

1.37 

0.062 

0.063 

16.5 

0.150 

0.50 

1.57 

0.047 

0.078 

18.9 

0.100 

0.60 

1.90 

0.032 

0.093 

22.8 

For  millimeter-wave  operation  in  the  vicinity  of  100  GHz,  the  following 


tables  provide  results  of  interest.  For  T  +  H  ^  0.002  inch 


Table  VIII.  (R  =  0.004  inch) 


(T/2R) 

(H/R) 

(27lR/\  ) 

O 

T(  IN) 

H(  IN) 

f  (CH  ) 

O  2 

0.250 

0. 

1.57 

0.00200 

0. 

738. 

0.225 

0.05 

1.59 

0.00180 

0.0002 

747. 

0.200 

0.  10 

1.60 

0.00160 

0.0004 

752. 

0.150 

0.20 

1.73 

0.00120 

0.0008 

813. 

Table  IX. 


(T/2R) 

(H/R) 

(27TR/A  ) 

0 

T(IN) 

H(IN) 

T+H(IN) 

R(IN) 

f  (GH  ) 

0  z 

0.250 

0. 

1.57 

0.014760 

0. 

0.014760 

0.02952 

100. 

0.225 

0.05 

1.58 

0.013446 

0.001494 

0.014940 

0.02988 

99.4 

0.200 

0.10 

1.60 

0.012032 

0.003008 

0.015140 

0.03008 

100. 

0.150 

0.20 

1.73 

0.009756 

0.006504 

0.016260 

0.03252 

100. 

0.113 

0.12 

2.23 

0.009469 

0.005028 

0.014497 

0.04190 

100. 

Table  X. 


(T/2R) 

(H/R) 

(2rcR/X  ) 

O 

T(IN) 

H(  IN) 

T+H(IN) 

R(IN) 

f  (GH  ) 

0  z 

0.025 

0.0 

6.9 

0.00648 

0. 

0.00648 

0.1296 

100. 

0.010 

0.0 

13.5 

0.00508 

0. 

0.00508 

0.2540 

99.9 

0.005 

0.0 

21.4 

0.00402 

0. 

0.00402 

0.4020 

100. 

0.00167 

0.0 

44.5 

0.00278 

0. 

0.00278 

0.8340 

100. 
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4 .  Discussion  and  Conclusions 

The  data  plotted  in  Figs.  4  and  5  were  generated  to  simulate  dielectric 
resonators  formed  in  GaAs  substrates.  For  these  data  the  height  of  the  upper 
ground  plane  above  the  resonator  structure  (see  Figs,  la  and  lb)  was  taken  to 
be  sufficiently  large  so  that  the  influence  of  the  ground  plane  on  the  circuit 
performance  would  be  insignificant.  This  height  was  determined  by  evaluating 
performance  for  several  different  height  values  for  each  resonant  structure  to 
determine  that  height  value  for  which  the  interaction  becomes  insignificant 
with  regard  to  influencing  the  first  three  significant  figures  in  computed 
values  of  T/2R. 

The  data  plotted  in  these  figures  is  applicable  to  resonators  operating 
at  frequencies  occurring  from  the  microwave  through  the  millimeter-wave 
regime.  It  should  be  noted  that,  although  the  analyses  described  in  this 
report  were  developed  to  be  invoked  for  groove-type  resonators  (with  wide 
groove  approximation),  they  are  applicable  to  and  even  more  accurate  for  all 
GaAs  resonator/substrate  combinations  where  a  circular-cylindrical  disk 
resides  on  a  substrate  (i.e.  no  grooves).  In  addition  to  the  information 
contained  in  Figs.  4  and  5,  Tables  I  through  VII  can  be  useful  for  microwave 
designs  of  such  structures.  Tables  VIII  through  X  are  similarly  useful  for 
millimeter-wave  realizations. 

Once  the  analyses  described  in  section  2  of  this  report  were  operational 
the  investigation  of  groove-type  resonators  became  one  of  satisfying  the 
following  conditions.  In  order  to  accurately  characterize  resonators 
including  coupling  effects  associated  with  the  groove  and  microstrip  feed  line 
interactions,  reasonable  (not  too  large  or  too  small)  geometries  needed  to  be 
determined  to  enable  experiments  to  be  performed  at  microwave  frquencies  where 
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automated  equipment  could  enable  accurate  comprehensive  resonator  characteri¬ 
zation.  Furthermore,  fabrication  of  versatile  test  samples  is  simplified  at 
microwave  frequencies  due  to  the  larger  physical  size  required,  thereby 
enabling  less  expensive  empirical  work  to  be  conducted.  On  the  other  hand, 
when  models  are  scaled  for  operation  at  95  to  100  GHz,  the  geometries  must  be 
compatible  with  processes  employed  in  monolithic  circuit  work. 

Based  upon  the  data  provided  in  Tables  I  through  X,  the  following  obser¬ 
vations  were  made.  In  order  to  fabricate  scaled  (microwave)  model  resonators 
in  K=13  material  having  resonant  frequencies  below  2  GHz,  excessively  large 
diameter  resonators  are  needed  (diameters  in  excess  of  10  inches).  Further¬ 
more,  by  virtue  of  the  resonant  frequency  being  affected  by  both  resonator 
diameter  and  thickness,  it  is  observed  that  for  100  GHz  geometries  (Tables  IX 
and  X)  that  either:  excessively  large  resonator  diameters  are  needed  for  GaAs 
substrates  which  are  then  enough  to  make  for  practical  reactive  ion  etching  of 
annular  grooves  (e.g.  0.168  inch  diameter,  0.0056  inch  thick  substrate;  or 
excessively  thick  substrates  for  practical  reactive  ion  etching  of  grooves 
with  reasonable  rsonator  diameters  (e.g.  0.0149  inch  thickness,  0.030  inch 
diameter).  Additionally,  for  100  GHz  resonators  the  coupling  microstrip  line 
used  to  launch  energy  to  the  resonator  is  estimated  to  require  a  GaAs  sub¬ 
strate  thickness  of  no  more  than  0.0027  to  0.005  inch  (see  Appendix  C)  in 
order  to  remain  free  of  higher  order  mode  propagation.  For  all  of  these 
reasons  it  is  felt  that  other  alternatives  to  dielectric  resonator  structures 
with  potential  compatibility  for  monolithic  realization  should  be  pursued  at 
this  point,  instead  of  the  groove-type  resonator  described  here.  Some  alter¬ 
natives  will  be  described  in  section  5.  It  should  be  noted  that  since  the 
alternatives  to  be  considered  in  section  5  do  not  require  groove-type 
processing,  the  microwave  experiments  needed  to  empirically  assess  the  effects 


of  the  resonator  groove  are  no  longer  needed.  Appendix  C  describes  some 
preliminary  work  performed  in  support  of  the  previously  anticipated  scaled 
experimental  efforts.  Hence,  the  effects  attributable  to  the  "wide  groove" 
approximation  are  no  longer  needed. 

5 .  Suggestions  for  Future  Work 

For  the  reasons  discussed  in  the  previous  section,  two  categories  of 
dielectric  resonators  should  be  considered.  One  category  encompasses  the 
formation  of  resonators  on  GaAs  substrates  by  processing  layers  of  higher 
dielectric  constant  materal  (e.g.  Barium  titanates)  on  the  substrate  surface. 
The  higher  dielectric  constant  material  should  reduce  the  geometries  needed  to 
achieve  mm-wave  resonance.  The  analyses  described  in  section  2  of  this  report 
are  flexible  enough  to  directly  treat  such  structures.  However,  significant 
effort  is  necessary  to  study  the  processing  methods  needed  for  realizations. 
A  variation  of  this  approach  is  to  sandwich  a  lower  dielectric  constant  layer 
(e.g.  k=2  to  4)  between  the  Barium  titanate  and  substrate.  This  could  raise 
the  Q. 

A  second  category  of  dielectric  resonator  which  might  warrant  further 
consideration  is  the  processing  of  K=13  material  on  GaAs  substrates.  The 
processing  methods  and  material  choices  need  further  consideration  and 
analysis.  This  category  is  not  currently  considered  as  promising  as  the  type 
mentioned  in  the  preceding  paragraph  because  the  geometry  reduction  is  not 
inherent  in  these  structures.  A  potential  advantage  might  be  accrued,  how¬ 
ever,  in  eliminating  the  need  for  dry  etching  thick  materials. 


APPENDIX  A  -  Program  DRES 


Program  DRES,  which  computes  the  characteristics  of  dielectric  circular 
cylindrical  resonators  (Figs,  la  and  lb)  without  fringing  is  written  in 
FORTRAN  IV  and  is  compatible  with  the  CDC  CYBER  175  digital  computer  system. 
With  minor  modifications  this  program  should  also  be  suitable  for  use  with 
digital  computer  accepting  FORTRAN  IV  encoded  programs.  Although  a  precise 
determination  of  the  execution  time  of  programs  has  not  been  determined, 
experience  has  shown  that  a  few  seconds  are  needed  to  determine  enough 
geometric  parameters  for  several  frequency  points  to  enable  plots  of  the  type 
shown  in  Figs.  2  through  5  to  be  made. 

This  program  embodies  the  analyses  described  in  section  2.1.  When  the 
proper  input  paramters  are  specific  this  program  evaluates  the  ratio  T/2R  (see 
Figs.  2  through  5)  for  given  values  of  27iR/\Q,  H/R,  D/R,  substrate  dielectric 
constant,  and  resonator  dielectric  constant.  The  input  parameters  required 
are  described  as  follows: 

ER  =  relative  permittivity  of  the  resonator  dielectric  material 
ES  =  relative  permittivity  of  the  substrate  dielectric  material 
FI  =  lowest  frequency  at  which  computations  are  to  be  performed  in  GHz 
F2  =  highest  frequency  at  which  computations  are  to  be  performed  in  GHz 
NK  =  number  of  frequency  intervals  between  FI  and  F2  to  be  used  in  compu¬ 
tations 

R  =  radius  of  dielectric  resonator  to  be  evaluated  in  inches 


The  output  parameters  provided  by  this  program  are  described  as  follows: 
A  =  evaluated  ratio  T/2R 
TWOPDL  =  ratio  2nR/k 

o 

D  =  ratio  H/R 
F  =  ratio  D/R 


ratio  D/R 
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The  following  should  be  taken  into  consideration  if  modifications  to  the 
program  coding  are  desired.  Lines  number  00310,  00320,  and  00330  are  used  to 
set  the  initial  value  (-0.024)  of  the  height  "H,"  the  number  of  steps  (11)  to 
be  taken  for  incrementing  the  initial  value  of  "H,"  and  the  size  of  the  step 
(0.025)  for  "H."  Mote  that  in  line  number  00330  the  first  computed  value  for 
D  is  0.001.  This  is  to  circumvent  problems  of  singular  behavior  incurred  for 
D=0 . 
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The  listing  for  program  DRES  follows. 


oo  i  o o  ppogprh  dfes  ■  i  nput  .  output  > 

001 10  PPINT  1 

00120  1 FOPNRT  ( 1  H  .♦EP.ES-F1.F2.NK  .PV) 

00150  PERB^.EP.ES.Fl  •  F  £  <  NH  .  P 

00151  Fl=Fl^l.E«-09 

00152  F£=F2*1 .E+09 

00 153  P=2.54E-0£vp 

00160  CSPEEB=3. OE+OR 

00170  FACT1  =  <2.  v3.  14 159>  -'C SPEED 

00 ISO  FRCT2P=  C£. 4 05*2 . 405>  s tR*P) 

00190  FRCT12=FRCT1*FRCT1 
00200  DELF=  <F£-F 1 > / 1 0 . 

00210  FPEQ1=F1-DELF 

00220  00  £  1=1.11 

00230  FPEG=FPEO I +LELF 

00240  FPE02=FPEQ^FPE*3 

0025  0  BRPG=FRCT  1  £vFPE€»2^EP-FPCT2P 

00251  BRPG=RBS<BRP6> 

0026  0  B=SOPT  <BRPG> 

00261  B=PvB 

00270  CRPG=FRCT2P-FRCT1£^FPEC*£ 

00271  CRPG=RBS  <CRPG> 

00230  C=SOPT <CRP6> 

00291  C=P*C 

00290  ERPG=FRCT2P-FRCT12*FPEQ£vES 
00291  ERPG=RBS  ■ (ERPG!> 

00300  E=SOPT<-'ERPG> 

00301  E=P>E 
00310  D I =-0.024 
00320  00  3  J= 1.11 
00330  D»D 1+0.025 
00331  D=Dv2.54E-02'P 
00340  FI =-0.2 49 
0035  0  DO  4  K  =  1 .  NK 
00360  F=FI+0.25 
00-361  F=F*2 . 54E— QSs'P 
00370  CD=OD 
00380  EF=E*F 
00390  HTRN1*TRNH.:CD> 

00400  HCTN1=1 ./HTRN1 
00410  HTfiN£=TPNH'EF> 

00420  HCTN£=1 . /HTRN2 
00430  RTRP61*  <C/B>  ♦HCTN1 
00440  RTRPG2=  >'E '*B>  *«CTN£ 

00450  R=  <\  .  /  <'2.  4B>  )  ♦<RTflN':flTfiP61>  +RTRN  'RTRPG£>  > 

0046 0  TWOPDLsFpi:  T 1  ♦FREOP 
00470  PPINf  5«R. TWOPDL.D.F 

00480  5FDPNRT ( 1H  . 2HR=.F7. 3. 3X. 7HTMDPDL=. F7. 3.3X» £HH=.F7. 3*  3X. 2HD=. F7. 3> 

00490  FI=FvPy 2. 54E-02 

00500  4CDNTINUE 

00510  DI*D4P^2.54E-02 

00520  3C0NTINUE 

00530  FPEQI =FPEO 

00540  2C0NTINUE 

00550  STOP 

00560  END 

PERDY. 
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APPENDIX  B  -  Program  BETACI 


Program  BETACI  is  written  to  determine  the  value  of  P  . ,  as  described  in 
section  2.2,  so  that  it  can  be  used  in  program  DRES  for  taking  into  account 
fringing  field  effects  of  the  dielectric  resonator  behavior.  Accessed  by 
program  BETACI,  but  not  listed  here,  are  four  subroutines  available  from  the 
Control  Data  Corporation,  Cybernet  Services,  International  Mathematical  and 
Statistical  Library  (IMSL).  These  subroutines  are  named:  MMBSJO,  MMBSJ1 , 
MMBSKO,  and  MMBSK1.  MMBSJO  computes  values  of  the  Bessel  function  of  the 
first  kind  of  order  zero.  MMBSJ1  computes  values  of  the  Bessel  function  of 
the  first  kind  order  one.  MMBSKO  computes  values  of  the  modified  Bessel 
function  of  the  second  kind  of  order  zero.  MMBSK1  computes  values  of  the 
modified  Bessel  function  of  the  second  kind  of  order  one. 

The  input  parameters  for  this  program  are  described  as  follows: 

FREQ  =  frequency  to  be  used  in  evaluation  in  hertz 
ER  =  relative  permittivity  of  the  dielectric  resonator  material 
BETA  =  initial  estime  to  be  used  in  search  routine  to  determine  Pci  (see 
explanation  to  follow) 

R  =  radius  of  resonator  in  inches 

To  determine  a  value  of  BETA  to  serve  as  an  input  parameter  consider  the 
following.  Since  the  argument  needed  for  square  root  computations  (e.g.  line 
number  00210)  must  be  non-negative,  then 

R2(— )2  (e  -l)f2  -  (P  .R)2  £  0 
c  r  rci 

This  means  that 


R2 


> 


(p  .R)2 

Kci 

(— )2  (e  -l)f2 
c  r 
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be  defined  by 


Let  R 


3.82 

(R 


(3 . 82)2 


3.82; 


(^?)2(er'1)f2 


Let  the  input  parameter  BETA  be  determined  by 
2.405 


BETA  = 


R 


3.82 


This  assures  that  the  argument  for  certain  roots  in  program  BETACI  are  posi¬ 
tive  with  occurring  for  the  range  from  2.405  to  3.82  (the  first  non¬ 
trivial  roots  for  JQ  and  3^ ,  respectively. 

Upon  receiving  the  input  parameters  program  BETACI  employs  a  Newton- 
Raphson  method  to  converge  to  a  value  of  0  with  accuracy  to  four  significant 
figures  based  on  the  convergence  criterion  used  in  line  number  00340. 

The  output  quantities  provided  by  program  BETACI  are  described  as  fol¬ 


lows  : 

BETAR  =  last  value  of  produce  PciR  determined  in  the  search  process 

BETAPR  =  the  last  value  of  PciR  determined  in  the  search  process  but 
increased  by  1/100%  for  use  in  finite  difference  approximations 
employed  in  convergence  procedure 

IERJ0  =  output  error  parameter  determined  in  routine  MMBSJO  (see  IMSL 
writeup  on  MMBSJO) 

IERJ1  =  output  error  parameter  determined  in  routine  MMBSJ1  (see  IMSL 
writeup  on  MMBSJ1) 

IERK0  =  output  error  parameter  determined  in  routine  MMBSK0  (see  IMSL 
writeup  on  MMBSK0) 

IERK1  =  output  error  parameter  determined  in  routine  MMBSK1  (see  IMSL 
writeup  on  MMBSK1) 

XNUM  =  numerator  used  in  evaluating  increment  of  for  Newton-Raphson 

method 
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XDEN  =  intermediate  quantity  used  for  debugging  (currently  super¬ 
fluous) 

XDENOM  =  denominator  used  in  evaluating  increment  of  8  ^R  for  Newton- 
Raphson  method 

It  should  be  noted  that  all  of  the  above  output  parameters  are  used  for 
observing  convergence  behavior  and  are  typically  not  necessary  as  output 
quantities  for  normal  convergence.  IERJO,  IERJ1,  IERXO,  and  IERK1  are  equal 
to  zero  for  properly  convergent  operation. 

Finally , 

BETCAL  =  calculated  end  value  of  needed  for  use  with  program  DRES  in 
(meters)  * 

Once  BETCAL  is  determined  from  program  BETACI,  this  value  is  then  used  in 
program  DRES  for  the  same  parameters  by  placing  line  number  00180  in  program 
DRES  with  "00180  FACT2R  =  BETCAL* BETCAL . "  Here,  the  actual  value  of  BETCAL 
would  be  used  in  lieu  of  the  name  BETCAL. 


The  listing  for  program  BETACI  follows. 


00 1  0  0 
0  o  to; 
00 1  05 
00110 
00 1  3  0 
00130 
00  14  0 
00141 
0015  0 

ooiro 
00131 
00133 
00133 
00134 
00135 
00130 
00131 
00300 
00301 
0031  0 
0031 1 
0033  0 
00331 
00330 
00331 
0034  0 
00341 
00350 
00351 
00360 
0037  0 
00330 
00390 
00  300 
00  301 
0031  0 
00315 
00330 
00331 
0033  0 
00331 
00333 
00340 
00350 
00  36  0 
PERCY 


PP06PPH  BETPCI -INPUT. OUTPUT} 

PERL  MM B'i  J 0 ,  MM  FT  J  t  . MMF'  K  0 •  MMB3K  1 
CPLL  IM3L 
PPINT  1 

1FOPMPT ( 1H  . *FPE9,EP,B£TP,P*> 

PEPO*.FPEO« EP. BETP.P 
P=P*3.54E-03 
BETPP=PETP*P 
I  OPT  = 1 

FPCT1  =  ('3.  *3.  14159  -3.  E+08>  *'-3.  *3.  14159-  3.  E*08> 

BETPPP=1. 0  0  0 1 *BETRP 
3 IF 'BETPPP. 6E. 3.33>4.5 
4PPIMT  6 

6F0PMPT  < 1 H  « ♦BETPP  PNC  BETPPP  TOO  LPPGE*> 

60  VO  3 

5X.J0*MMBi  JO  <BETPP.  IEPJ0> 

XJ0P*MMB3J0  1  BETPPP-.  IEPJ0> 

X J 1 =MMBSJ 1 -BETPP. I EP  J 1 > 

X J1P-HHBSJ1 'BETPPP. I EP  J 1 1 

UPPG*  ( P  *P*FPCT 1 ♦  •  EP- 1 .  ■  *FPEO*FPE©:>  -  ':BETPP*BETPP> 

UPP6P*  .:P*P*FPCT1  *•  EP-1 .  >  *FPEO*FPEQ>  -< BETBPP*BETPPP> 
U=S0RT<UBP6> 

UP=3QPT  (IJPP6P) 

XK0=MMB  iK  0  •:  IOPT.  U.  IEPKO'j- 
XK  0P*MMBSK  0  -'I  OPT  .UP.  IEPK  0> 

XK 1 =MMP  SK 1 < I OP  T, U , I EPK 1 > 

xkip=mmb$ki  <iopt,up»  iepkd 

XNUM*  ■ BETPP* 'XJO-'XJt  >  >  ♦  dj* •'  XK  O.XK  1>  > 

XMUMP* '  BETPPP*  ■' X JOP' X  J 1  P>  »  *  <UP* e XK OPxXK  IP}  > 

T 1  =X.J  0-X.J  1 

T3*FETPP*c  f  ''X  JO-'XJl  '  ♦  1 . /BETPP}  -  •:  < X J OxXJ  1  .>  **3 .  )  }• 

T3»3 .  ♦  ■  BETPP  -1.0  ♦  ('XK  0  -XK' 1  > 

T4=BETPP*  < (XK 0  - XK 1 ) **3 . > 

XDEN=T1 +T2-T3+T4 

XDENOM*  < X M IJ M P — X N U M >  --  ('EETPPP-BETPR} 

H=-XNUM.-XDENOM 
PH*BBS (H) 

BETBP=BETBP+H 
BETPPP* 1. 0001*BETBP 

PPINT*. BETPP. BETPPP, IEPJ0.IEPJ1* IEPKO. IEPK 1 , XNUM. XBEN. XDENOM 

BETCPt.=  BETPP-P 

PPINT*. BETCPL 

I F ( PM . L  T . l.E-04} 3*3 

3  STOP 

END 
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APPENDIX  C  -  Preliminary  Work  for  Experiments 


The  data  listed  in  Tables  I  through  VII  were  evaluated  to  provide  infor¬ 
mation  needed  for  guiding  geometries  to  be  used  in  experiments  frequency 
scaled  to  microwave  frequencies  where  automated  equipment  could  enable 
accurate,  comprehensive  resonator  characterization.  It  was  initially  envi¬ 
sioned  that  experimental  realizations  such  as  those  shown  in  Figs.  6,  7,  and  8 
could  be  used.  Figure  6  shows  a  "groove-type"  resonator  fabricated  in  TiO^- 
loaded  stycast  with  a  microstrip  line  coupled  to  the  resonator.  The  rf  signal 
is  launched  and  measured  through  N-type  coaxial  connectors  as  shown  in  Fig.  6. 
In  Fig.  6  the  top  ground  plane  is  removed  to  expose  the  circuit  nature. 
Figure  7  shows  the  two  parts  used  to  form  the  resonator-circuit  structures. 
The  plug  seen  in  this  figure  is  machined  as  well  as  the  counter-bored  hole  in 
the  substrate,  into  which  the  plug  is  set  and  affixed  with  Emerson  and  Cuming 
stycast  loaded  epoxy  having  a  dielectric  constant  of  13  to  match  the  dielec¬ 
tric  constant  of  the  stycast  board.  The  microstrip  line  and  transitions  to 
N-type  coaxial  connectors  were  optimized  for  match  on  these  boards  using  a 
Hewlett  Packard  1815A  time-domain  ref lectometer .  The  assembly  ultimately  used 
in  early  experiments  employed  an  upper  ground  plane  and  is  shown  in  Fig.  8. 

Preliminary  measurements  were  needed  on  isolated  microstrip  lines  fabri¬ 
cated  on  stycast  boards  to  ensure  that  connectors,  transitions,  and  microstrip 
lines  were  functioning  adequately  up  to  2  or  3  GHz.  Initial  measurements  were 
performed  using  stycast  boards  (k=13)  having  a  thickness  of  *5  inch.  However, 
measurements  showed  that  reasonable  transition  performance  for  straight 
through  microstrip  lines  could  not  be  maintained  above  about  1.2  to  1.5  GHz. 
Upon  making  similar  measurements  using  h,  inch  thick  stycast  boards  it  was 
found  that  adequate  performance  could  be  achieved  to  above  2  GHz.  Figures  9 
and  10  show  the  insertion  loss  and  VSWR  measurements  performed  on  the  k  inch 


thick  board.  The  limitation  on  substrate  thickness  observed  by  comparing 
*5  inch  thick  to  \  inch  thick  board  results  is  attributed  to  that  frequency  at 
which  higher  order  modes  are  launched  in  the  microstrip  line. 

Prior  to  that  point  in  time  that  accurate  analyses  were  ready  for  use, 
gross  approximations  were  used  to  estimate  geometries  for  experimental  models 
so  that  coupling  effects  could  be  studied.  It  was  found  that  appropriate 
resonator  geometries  were  not  well  estimated  which  resulted  in  resonators  not 
having  detectable  resonant  frequencies  in  the  frequency  range  capable  of  being 
measured.  Once  the  analyses  described  in  section  2  of  this  report  were  opera¬ 
tional  the  problems  encountered  in  early  experiments  were  better  understood. 


Fig.  6  —  Photograph  of  groove-type  dielectric  resonator  coupled  to  microstrip  line  fabricated  in  k  =  13  Stycast  material. 
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Fig.  9  —  Measured  insertion  loss  versus  frequency  for  back-to-back  N-type 
to  microstrip  (nominally  50-ohm)  transitions  on  1/4  inch  thick  k  *  13 
Sty  cast  board. 
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Fig.  10  —  Measured  VSWR  versus  frequency  for  back-to-back  N-type  to 
microstrip  (nominally  50-ohm)  transitions  on  1/4  inch  thick  k  -  13  Stycast 
board. 
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